Rising time of entanglement between scattering spins 
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We investigate the time evolution of entanglement in a process where a mobile particle is scattered 
by static spins. We show that entanglement increases monotonically during a transient and then 
saturates to a steady-state value. For a quasi-monochromatic mobile particle, the transient time 
depends only on the group-velocity and width of the incoming wavepacket and is insensitive to the 
interaction strength and spin-number of the scattering particles. These features do not depend on 
the interaction model and can be seen in various physical settings. 
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Scattering, an almost ubiquitous mechanism in physics 
and a rather well-studied topic, has very recently gained 
renewed interest within the community working on quan- 
tum mechanics in virtue of the potential that it has in 
many respects^* 2 -'^. Scattering processes between two 
subsystems are effective in order to probe correlation 
properties of quantum many-body systems^. Moreover, 
under proper conditions, scattering can be exploited to 
prepare non-classical states of unaccessible systems 2 -^. 
Interesting proposals have been put forward for the gen- 
eration of quantum correlated light-matter states via off- 
resonant coherent scattering 5 . Frequently, a stationary 
approach to scattering processes is chosen, especially for 
the sake of entanglement generation: the system under 
scrutiny is observed at long-time scales, when it should 
have reached steady conditions. This cuts the (often com- 
plicated) time-evolution from the effective description of 
scattering dynamics. While this approach is computa- 
tionally convenient and frequently useful, it is not fully 
satisfactory since it leaves some interesting questions, re- 
lated to the details of the dynamical evolution, unan- 
swered^. In particular, it does not give information on 
the time needed by entanglement to reach its stationary 
value, which is a pivotal point for the aims of coherent 
quantum information processing (QIP). In fact, an esti- 
mate of such transient time will help us to anticipate, 
quantify and eventually counteract the effects that deco- 
herence might have in a given process. 

In this paper, we focus on a prototypal setting involv- 
ing both mobile and static spin particles 2 ^^ to shed some 
light onto these issues. We clearly identify the physical 
parameters that determine the duration of the scattering 
process. Independently of the Hamiltonian chosen in or- 
der to model the spin-spin coupling, we show that when 
a quasi-monochromatic incoming wavepacket of the mo- 
bile spin is prepared, with a given average momentum 
feo, stationary conditions are reached in a time dictated 
only by the wavepacket width in momentum space Ak. 
This parameter can thus be used to tune the duration of 
the scattering event and the rate of entanglement gener- 
ation in the system. Counter-intuitively, the interaction 



strength of a given spin-spin Hamiltonian model coupling 
mobile and static particles does not affect the scattering 
transient time At but only determines the maximum en- 
tanglement attainable in the process. Although, as an- 
ticipated, our conclusions do not depend on the details 
of the model being considered, most of our quantitative 
results are presented for spin-spin interactions induced 
by the Heisenberg exchange coupling. This model arises 
naturally in the spin-spin interaction of magnetic impuri- 
ties embedded in a one-dimensional electron nanowire as 
well as in other situations, including the interaction be- 
tween a free and a bound electron in a carbon nanotube 2 . 
However, in the second part of our work we quantitatively 
address the results for an anisotropic XYZ model. Tun- 
ing the ratio of its parameters, we are able to span a wide 
range of significant spin-spin Hamiltonians. Remarkably, 
our results can be extended to an arbitrary number of 
static particles and intrinsic spin numbers, which makes 
them valid in a heterogeneous set of physical settings 
(ranging from spintronics to cavity-quantum electrody- 
namics )2£d&l. In passing, we also reveal an unexpected 
and general monotonic time-behavior of entanglement, 
observable under easily-matched conditions. Although 
wiggling of the reflected wavefunction of the mobile par- 
ticle is observed due to interference at the scattering 
site(s), entanglement can only grow in time. First, we 
empirically observe this behaviour for relevant coupling 
models and then give a clear physical explanation for the 
case of the Heisenberg exchange coupling and one scat- 
tering center. The features described above make it clear 
that a scattering-based mechanism holds the promises for 
a genuinely control-limited distribution of entanglement 
in a partially addressable quantum register. In contrast 
to procedures based on the temporal gating of spin-spin 
interactions the process addressed in this paper allows 
for a dramatic relaxation of the time control, which is 
a major advantage in light of a potential experimental 
implementation. 

The remainder of the paper is organized as follows. In 
Sec. U we introduce the single-scatterer version of the sys- 
tem addressed in our work and study the entanglement 
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rising-time. An important bench-mark is set by com- 
paring the results of this study with those corresponding 
to the scattering-less case of two exchange-coupled static 
particles. Sec.HTlextends the analysis to the two-scatterer 
case. In Sec. IHII we quantitatively prove the insensitivity 
of the entanglement rising-time to the spin-spins coupling 
strength. Finally, in Sec. IIVI we summarize our findings 
and remark their most important implications. 



I. SCATTERING BY A SINGLE STATIC 
PARTICLE 

In order to provide a significant mile-stone for our main 
results, we first review the well-known case of two static 
spin-1/2 particles, e and 1, interacting via the Heisenberg 
coupling H = J&-Si. Here, a (Si) is the spin operator 
of particle e (1) and J is the coupling rate (we use units 
such that h = 1 throughout the paper) . We assume the 
initial state 

lx(0)) = |U) el = ^(|* + ) + |*-))el, (1) 

where |T), \i) are the spin states of each particle and 
|1r±) el = (|T|> ± ||T))ei/V2 are steady states of H with 
energies J /4 and —3^/4, respectively. Upon evolution, 
the e— 1 state reads 

| X (r)> = ^ ljT,i \^> + )+e ajT ' i |*-» e i. (2) 

As a measure of the entanglement between the particles 
in Pei(T~) = \x{t))(x(t)\, we use the logarithmic negativ- 
ity Ejs[(t)£. A straightforward calculation gives En(t) = 
log 2 (1+| sin Jt\ ) , which oscillates with characteristic time 
J- 1 . E n (t) is maximized at r=(2g+l)7r/2 l 7- 1 (q € N), 
which shows that for a set coupling strength J a fine 
control on r is required. This is not the case when the 
spin-spin interaction takes place during scattering. 

To show it, let us address the case of a mobile parti- 
cle e, interacting with the static spin 1 during scatter- 
ing. To fix the ideas, we consider a one-dimensional wire 
along the ir-axis where spin 1 is embedded at position 
x = 0. Spin-spin interaction occurs, via a Heisenberg 
model, when the mobile spin reaches x = and is scat- 
tered by 1. Assuming a quadratic dispersion law for e, 
the Hamiltonian reads 

H=^— + J&-S 1 S(x), (3) 
2m* 

where p ( m*) is the momentum operator (effective mass) 
of e and J is the coupling strength (notice that, unlike 
J, J has the dimensions of a frequency times length). 
The incoming mobile spin e is supposed to have wavevec- 
tor r\k (r) = ±1, k > 0) and be prepared in the spin 
state =!, t) e while the static spin 1 is initially in 
l m i = .UT)i- We call fi — {Tn' e7 m[} {v = {m e ,mi}) the 



set of spin quantum numbers of mobile and static parti- 
cles before (after) the scattering event. Correspondingly, 
r kriv (*fcni/) ^ S t ne probability amplitude that e is re- 
flected (transmitted) in state \m e =[, |) e , while the spin- 
state of 1 is changed into |toi ={, |) r These probability 
amplitudes depend on J/v^ (vk = k/m*) and are com- 
puted by solving the time-independent Schrodinger equa- 
tion (SE) associated with Eq. ([3]) and imposing proper 
boundary contitions at x — 0^-. The steady state of the 
system ) has wavefunction ^(x, v) — (x,v\ty£ v ) 
of the form 

(4) 

where we have omitted a factor 1/V27T and 9(x) is the 
Heaviside step function. As the kinetic energy of e is 
the only free-energy term in Eq. ^ the system's spec- 
trum is £fc = fc 2 /(2m*) and thus coincides with that of 
a free-propagating e. In order to investigate the dynam- 
ics of entanglement during scattering, we consider e as 
prepared in a Gaussian wavepacket \(p) such that 

(x\<p) = tp(x) = (2irp)- 1 / 4 e lko{x+Xo) e- ix+Xo)2/ ^ (5) 
with xq-, fco > the average position and momentum of 
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FIG. 1: (Color online) a) Em{t) (full curves) and free- 
propagating-particle / e (0, t) (dashed lines, rescaled to the re- 
spective maximum) in a log-lin scale for xo = 5 Ax, J/vk = 1 
and Afc/fco = 1(T 4 , 1(T 3 and 10~ 2 . b) E n (t) for x = 5Ax, 
Afc/fco = 10~ 2 and J/vk = 1,3, 10. The top line shows the 
same curves rescaled to their maximum. Fhe choice of the 
parameters used in these plots optimize their visibility (other 
values are equally valid). 
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e. As for the uncertainties, we have Ax — (Afc) 1 = v / 2/3- 
The overall initial state is taken as |^(0)) = \ ip) |T, |) e 1 , 
whose time-evolved version |^(t)) is found by solving the 
respective time-dependent SE. 

In Fig. [TJ we study the evolution of spin-spin entan- 
glement between e and 1 when the former is prepared in 
a quasi-monochromatic wavepacket. Clearly, in the pro- 
cesses at hand, the entanglement associated with the full 
quantum state involves both motional and spin degrees of 
freedom. The study of this hybrid form of entanglement, 
which is in general a difficult task requiring ad hoc tools 
for his approach, is beyond the scopes of our work. A pos- 
sible way to tackle it would be the investigation of the 
non-locality content of the state of particles e and 1 and 
involving both internal and external degrees of freedom, 
along the lines as in Rcf. 11 . In Fig.QJa) we set J/vk = 1, 
xq = 5Ax and plot En(t) for increasing values of Afc/fco 
up to 10 -2 . The differences with respect two static spins 
are striking. First, there is no oscillatory behavior. At 
small t, no entanglement is found as far as e has not 
yet reached 1. At large times, e is far from spin 1 and 
En{t) takes a steady value. At intermediate times, the 
entanglement shows a monotonic increase before eventu- 
ally saturating to a steady value, which does not depend 
on Afc. Second, unlike the case of static spins, the char- 
acteristic time At of the entanglement evolution is now 
independent of the spin-spin coupling strength J. To il- 
lustrate this, in Fig. [ljb) we set Afc and study En(t) 
for various J/vk 's. Remarkably, although the steady 
value of En(t) depends on J/vk , the rising time At does 
not, as it is revealed by the top-most curve in Fig. Hfb). 
There, for an assigned value of J/vk , we have rescaled 
each £jv(t) to the respective stationary value and found 
that the curves are identical. This marks a profound dif- 
ference with the static-particle case: the introduction of 
motional degrees of freedom in our problem of interact- 
ing spins does not result in a mere spoiling effect of the 
system spin coherences, as it might be naively expected, 
but deeply affects the dynamics of the particles involved 
in the process. To identify the parameters upon which 
At depends in the dynamical-scattering process, we first 
observe that in Fig. [TJa) an increase in Afc of one order 
of magnitude makes At ten times smaller, suggesting an 
inverse proportionality between these quantities. Then, 
in Fig.[lja) (dashed lines) we plot the probability density 
f e (x = 0, t) of finding particle e at x = for J = 0, i. e. the 
free-propagation case. It can be clearly seen that the time 
needed by En (t) to reach its steady value coincides with 
the time during which e is found at x = with non- 
negligible probability in the absence of any scattering. 
Thus, by using the free-particle time-energy uncertainty 
principle^ we obtain At ~ l/(vk Ak), which explains 13 
the aforementioned inverse proportionality between At 
and Afc for a given fco In other words, in the scattering 
case and for a quasi-monochromatic mobile particle the 
characteristic time over which entanglement changes is 
decided only by the kinetic parameters specifying the in- 
coming wavepacket ip(x). As for a given J/vk the steady 



value of Em{t) is insensitive to Afc [cfr. Fig. 1(a)] we 
conclude that the rising time At can be tuned simply 
by adjusting Afc, with no effect on the stationary value 
of entanglement (for a set J/ Vk ) ■ This result is key to 
the study of the feasibility of scattering-based QIP pro- 
tocols^ 3 -^ and the quest for effective ways to relax the 
control on a system. We have checked that these fea- 
tures hold even when particle 1 has spin number s> 1/2. 
We point out that although Fig. 1 might at first glance 
suggest that entanglement rises linearly with time in gen- 
eral this is not the case. Rather, our findings show quite 
clearly that such functional behavior is dictated by the 
shape of the incoming wavepacket of particle e (here as- 
sumed to be Gaussian). 

II. SCATTERING BY TWO STATIC 
PARTICLES 

We now consider the situation where the mobile spin 
e is scattered by two static spin-1/2 particles, 1 and 2, 
placed at x = and x = d, respectively. The Hamil- 
tonian is the same as in Eq. ([3]) with the inclusion of 
the term Ja S2 S(x — d). Each steady state (// = 

{m^m^m^}) has energy £fc = fc 2 /(2m*) and differs from 
Eq. dH) for the replacement 6(r)x)*->9[r)x— (r)+l)d/2\ and 
the addition of (A% v u e ir > kx +B^ v e~ ir i kx )[9{x)-8(x-d)}. 
The coefficients 7^ ?; v (^y = A, B, r, t) depend implicitly on 
J/vk and kd. As shown ir* 3 -^, once boundary conditions 
at x = 0, d are imposed, the steady states can be deter- 
mined. This configuration has recently been proposed 
as a potential way to set entanglement between remote 
spins^ 3 -^ via mediation of e and only mild time-control: 
It is enough to wait for a time r>Ar, so that a steady 
state is reached. Clearly, the quantification of At in 
the case of an incoming wavepacket (instead of a plane 
wave as in Refs^i 3 -^) is key to estimate the influence of 
dissipation and decoherence affecting the scattering par- 
ticles. This would be a pivotal point in view of potential 
experimental implementations. We take the initial state 
|*(0)) = \(p,t) e ||,|) 12 , whose time-evolution |*(t)) is 
computed through the time-dependent SE. We focus on 
the dynamics of entanglement -E/v[/?i2(t)] between static 
particles 1 and 2. Here 

Pu{t) = Y^ / dx{x 1 m, e \^{T)){^{T)\x 1 m e ) (6) 

is the 1—2 state obtained by tracing over the spatial and 
spin degrees of freedom of e. We call 

p e (Q,r)= / dxf e (x,r) (7) 
Jo 

the probability to find e within the interaction region 
O = {x : < x < d} at time t. In Fig. [DJa) we set 
xq = 5 Ax, J — Vk , kod — it and study the behavior 
of En[pi2(t)] and /ree-propagating-particle p e (ri, t) for 
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FIG. 2: (Color online) a) E n [ Pi2 (t)] for Ak/k = 1(T 2 , 1CT 3 
and 10~ 4 for J/vk = 1 and kod = ty (solid lines) and free- 
propagating-particle p e (Zl,r) (dashed lines, rescaled to their 
maximum), b) £W[pi2(i)] for Afc/fco = 10 -2 , kod — ty, 
J/vk = 0.5, 1 and 3. The choice of the coupling strengths 
used in these plots optimize their visibility (other values are 
equally valid). 



increasing values of Afc/fco up to 10~ 2 . On the other 
hand, Fig. [2Kb) considers the effect of different choices of 
J/vk over the entanglement dynamics, for set values of 
Afc/fco and kgd. Features similar to those revealed for a 
single scattering center are found [see Figs.[TJa) and (b)]. 
Entanglement increases with time to reach a steady value 
that depends on J/vk and kod. As in the single-scatterer 
case, At is determined only by the time needed by the 
free-propagating (p(x) to cross 57. In Fig. [2] our focusing 
on the regime kod = ir is due to the fact that such setting 
allows for efficient entanglement distribution schemes 3 ' 4 . 
However, we have checked that other values lead to sim- 
ilar conclusions. 

In reaching our conclusions, we have performed a few 
technical steps which is instructive addressing here. In 
each of the studied configurations, the time evolution of 
1^(0)) = \<p) |/2) (with 1/2) the total initial spin stated) can 
be expanded over the basis of steady states )} with 

coefficients given by |\I>(0)). Provided that xq > 

3 Ax, i.e. at r = spin e is out of the domain x>0, these 
scalar products are well-approximated by taking 6{—x) ~ 
1 and neglecting contributions from other ^-functions in 



Eq. (JU), so that 

(% v mo)) CpM-*) 

with <p(k') the Fourier transform of (p(x) (an analogous 
expression holds in the many-scatterer case). Eq. ([5]) 
allows us to evaluate 

/>OC 

(x,»\*(t))=J2 / *(*^|*(0)) e - ie ^^ ;7) (x, M ), 

(9) 

which can be solved analytically via a power-series ex- 
pansion of each spin-dependent amplitude 7^ v (7 = 
A,B,r,t) around the carrier wavevector fco. The re- 
placement of the ro-th order expansion of these coef- 
ficients in (x, /J, \&{t)) results in integrals of the form 

J °° dke- ak2 ^ k ~k m (m = 0,.. > n,k = k/k and Re a >0), 
which can be computed in terms of exponential and er- 
ror functions. The behavior of f e (x,r) is then found by 
tracing the overall state p(r) = \^{r)) (&{t)\ over the spin 
degrees of freedom. As expected^, we find that during 
scattering f e (x,r) wiggles at x < due to interference 
between incoming and reflected waves. The spin state is 
found by tracing p(r) over the spatial and, in the many- 
scatterer case, spin degrees of freedom of ei£. 

III. MODEL-INDEPENDENCE OF 
ENTANGLEMENT RISING-TIME 

We are now in a position to explain the insensitivity of 
At to the values of J. Key to this task is the observation 
that as long as the wavepacket is narrow enough around 
fco the term proportional to <p(—k) (with fc>0) in Eq. (|S]) 
can be neglected. Thus C#£J*(0)) ~ 6 Vt+ 6p, tlx tp(k) 
meaning that, within the limits of our study, the spec- 
tral decomposition of 1^(0)) is the one corresponding to 
J = 0. On the other hand, the spectrum ek does not 
depend on J [see discussion after Eq. Q] so that the en- 
ergy and time uncertainties— Ae and At, are the same 
as in the free-particle case, i. e. At = 1/Ae ~ 1/ (vk Afc) . 
These features hold for a quasi-monochromatic incom- 
ing wavepacket. Most importantly, it is clear that our 
proof does not rely on the specific form of the interaction 
Hamiltonian, the number of scattering particles as well 
as their intrinsic spin numbers. While all these parame- 
ters affect the shape of r£ 's and t% ^'s, we have just 
shown that, for quasi-monochromatic wavepackets, they 
cannot influence the entanglement rising-time. 

To further illustrate the discussed insensitivity of At, 
in Fig. 3(a) we set Afc/fco = 1CP 2 and analyze the en- 
tanglement between 1 and 2 when these are coupled to 
e with strengths J e i=/= J e 2- The case of spin-1 scattering 
centers and equal couplings is also reported. In Fig. 3(b) 
we set the same Afc/fc as in Fig. 3(a) and address the 
XYZ spin-spin model 

£xyz= Jl[viSi,i6(x)+aiS 2tl 5(x-d)}. (10) 

l=x,y,z 
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We study the cases of J x = J y = J z /2 — Vk (embodying 
an XXZ model), J x = J y = Vk with J z — (isotropic 
XY) and J x = J y /2 — 3vk with J z — (anisotropic 
XY). The ratios of the parameters are chosen so as to 
provide the best visibility of each plots. The insets in 
Figs. 3(a) and (b) report all the curves rescaled to their 
maximum value, showing that the entanglement rising- 
time is not affected by the specific quantum spin-number 
or the interaction model, which only affect the station- 
ary value of entanglement. In particular, the applicabil- 
ity of our results to the XY model is remarkable since 
an effective XY model is indeed found^ considering the 
dispersive interaction of a single photon travelling across 
a ID waveguide (a GaAs/GaN nanowire, for instance) 
and static atom-like systems (such as InAs/GalnN quan- 
tum dots or nitrogen- vacancy centers in diamond^). The 
(typical) case of linear dispersion law for a photon cross- 
ing a waveguide matches the requirements of our study, 
i.e. a narrow-bandwidth mobile particle. These consid- 
erations demonstrate the broad applicability of our in- 
vestigation and results, which cover a wide range of ex- 
perimental situations, from spintronics to quantum op- 
tic o 2 '?' 4 '?'? . In the second scenario, in particular, the 
proven ability to experimentally engineer the tempo- 
ral shape of photonic wavepackets^ makes scattering- 
based techniques analogous to the one discussed here and 
in Refs^^ quite advantageous against strategies using 
time-dependent "modulation" or "pulsing" of spin-spin 
interactions. In the former case, once the kinematic prop- 
erties of the mobile-particle wavepackets are set, there is 
no necessity for time-control of the entanglement evo- 
lution. We remark that, in order for the steady-state 
entanglement not to depend on Ak [as in Figs. [1] and [2], 
this has to be smaller than the inverse of the character- 
istic length associated with 7-amplitudes. In particular, 
for two static spins, it must be Afc <l/d. 

For the Heisenberg model, the monotonic rise of entan- 
glement^ can be interpreted as due to the progressive 
construction of a phase-difference between spin compo- 
nents of the overall wavefunction. Here, we discuss the 
steps required in order to prove this feature. For the 
sake of argument, we focus on the case of Eq. and 
we assume that the incoming wavepacket ip(x) is quasi- 
monochromatic. The Hamiltonian H in Eq. ([3]) com- 
mutes with S 2 and S z , where S = <r + Si is the total 
spin of e and 1. It follows that for the initial spin state 
l*(0)) = |v)lT4) e ,i [see also Eq. ©] 



reflection and transmission probability amplitudes are 



1 



'el 



(11) 



where ip± (x, t) are evolved wavepackets fulfilling the con- 
dition <p a (x,0) =<p(x) and |*±) el = (|U)±||T)) e i/V2. 
Upon trace over the spatial variable of e we get En(t) — 
log 2 (l + |ImS(r)|) with £(t) the spatial overlap between 
<p~(x, t) and Lp + (x,r). Now, for |W ± ) el the spin-spin in- 
teraction in Eq. Q reduces to an effective static potential 
T±6(x) with T+ = J/4 and T_ = -3 J/4iS with associated 
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x/i+TVfT? 



D i&-[ccot(T a /v k ) 



(12) 



The steady states \^f v ) are similar to those in Eq. (|4|) 
with associated energies Sk — fc 2 /(2m*). Assuming 
x > 3Ax, the projections (^J^) of ip^fa) onto the 
stationary states are well-approximated by a form anal- 
ogous to Eq. (H]) (with due replacements). Provided that 
the width of the incoming wavepacket around fco is nar- 
row enough, we find (^ , fc 7) |<^ ± ) — $n,+ <f{k). Moreover, 



[t 



= r [*fc - *fco = * ]■ Moreover > 



± - *± 



t^]. Thus, by calling A r (A t ) the phase- 



difference between r + and r (t + and t ), we get 



E(r)^ / dx[\^ R (^r)\ 2 + \r+r-\e^\^ LL (x,r)\ 2 } 
+ / dx\t+t-\e iA *\ip ftR (x,T)\ 2 , 



(13) 

where ipf jP (x,T) (p — R,L) are evolved free wavepack- 
ets with average positions x p — s p xo at r — and 
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FIG. 3: (Color online) a) Two-static-spin case, Heinsenberg 
model. From top to bottom: En [pi 2 (t)] for J/vk ~ 1 (equal 
couplings, spin-1/2), unequal couplings J e i = 2J e 2 = 2.6«fc 
(spin- 1/2) and J/«k =1 for spin— 1 scattering centers (equal 
couplings), b) Two-static-spin case, spin-1/2, XYZ model. 
From top to bottom: J x — J y /2 = 3vk with J z = (anisotropic 



XY), 
J z ft = 



-v k . 



with J z = (isotropic XY) and J x = Jy 



(XXZ). Insets: curves rescaled to their maximum. 



We have taken Afc/fco = 10 and kod = w. 
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Fourier transforms tp(—s p k') with sl = — sr = 1. In 
our derivation we have neglected terms proportional 
to J dx tp*j R {x,T)ip f^(x,r) due to the assumption of 
dealing with an incoming wavepacket narrow enough 
around fco- As Lpf.L{x,r) and tpf^{x,T) are left- 
and right-propagating wavepackets, respectively, both 
J_ oo dx\ipf^L(x,T)\ 2 and f^°dx\(pf,R(x, t)| 2 are increasing 
functions of time whose shapes are similar to those 
in Fig. 2(b) for a Gaussian wavepacket <p(x). Hence, the 
behavior of |ImE(r)| depends on sinA r ( t ). As A r( - t ) S 
[0, 7r/2] (regardless of J and kg) the sine functions are al- 
ways positive. Therefore, |IulI](t)| is an increasing func- 
tion of time, proving the monotonic character of En(t). 
Physically, it is now clear that the build-up of entangle- 
ment relies on the phase-difference acquired by the singlet 
and triplet components of the incoming wavepacket once 
e is scattered off. Dynamically, as soon as the scattering 
process starts and the reflected and transmitted waves 
are progressively generated, the singlet and triplet com- 
ponents start to build-up the mentioned phase-difference 
and entanglement grows. For J/vk up to ~ 10 and 
Afc/fco up to 1CP 2 , the largest relative increments of 
and t k in the range [fc — 3A/c, fc +3Afc] are of the order of 
~ 10%, confirming the validity of our arguments. Clearly, 
the monotonic increase of entanglement does not hold for 
any initial spin state. For instance, starting from an en- 
tangled state of e and 1, Em may decrease with time. 

The study of wavepackets having Ax — ► 0, namely 
the regime opposed to quasi-monochromaticity, is under 
ongoing investigations. 



IV. CONCLUSIONS 

Our work reveals that in a scattering process involv- 
ing one mobile and many static spins, a situation that 



can be engineered in various physical settings2»2»£iiiii the 
time required for entanglement to reach its steady value 
can be tuned by preparing the mobile spin's wavepacket, 
regardless of the specific spin-spin interaction model. If 
a wavepacket narrow enough in frequency is prepared, 
the strength of the interaction between mobile and static 
spins determines just the stationary value of the entan- 
glement which, while scattering takes place, can only 
grow in time. Our findings show that the degree of 
control on the class of systems we have addressed can 
be significantly reduced down to the off-line engineering 
of a proper mobile spin's wavepacket. The quantifica- 
tion of the scattering transient-time enables the analy- 
sis of noise effects, such as phase-damping affecting the 
static spins, on protocols for entanglement distribution 
via scattering 2 -^. 
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